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Abstract 

We give an introduction to the techniques from microlocal analysis that 
have successfuhy been apphed in the investigation of Hadamard states of free 
quantum field theories on curved spacetimes. The calculation of the wave front 
set of the two point function of the free Klein-Gordon field in a Hadamard state 
is reviewed, and the polarization set of a Hadamard two point function of the 
free Dirac field on a curved spacetime is calculated. 
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1 Introduction 



The general framework of this work is that of quantum field theory on curved space- 
time. Since we still do not have a full theory of quantum gravity, as a first step we 
try to include gravitational effects into the quantum field theoretical description of 
our world by a semi-classical approximation. The gravitational field is considered as 
a fixed classical background spacetime on which matter is described by a quantum 
field theory. Since the Planck length, at which effects from quantum gravity are 
expected to become important, is very small, this approximation should have a wide 
range of validity and should be appropriate to predict quantum effects near black 
holes or in the early universe. One famous result in this framework is the Hawking 
radiation of black holes. 

The main problem in the quantization of a field theory on a curved background 
is that one cannot make use of Poincare invariance like on Minkowski spacetime, 
where this symmetry fixes a unique vacuum state. On a curved spacetime, however, 
there is no such natural choice of a state, and the central task is to single out the 
physically acceptable ones. 

For free, linear models it has turned out that these acceptable states are the 
quasi-free Hadamard states. They are determined by their two point function which 
has got a prescribed singularity structure fixed by geometry. For the free Klein- 
Gordon field it was recently discovered by Radzikowski that this Hadamard condition, 
which at first sight is a global condition, can be reformulated as a local condition 
in terms of wave front sets. This allows one to apply the elegant and powerful 
mathematical techniques from microlocal analysis. In this way it could be proved 
that certain states known for a long time are indeed Hadamard states. A microlocal 
spectrum condition for interacting quantum field theories was formulated and used 
to define Wick polynomials and to prove the renormalizability of interacting scalar 
field theories on globally hyperbolic spacetimes. 

For more realistic quantum field theories, however, first of all the Dirac field, 
more or less nothing was done, not because of conceptual problems, but because of 
the greater complexity due to the multi-component nature of the Dirac field. In this 
work the singularity structure of the two point function of the free Dirac field will 
be investigated. We will make use of the polarization set which was introduced by 
Dencker in order to analyze the singularities of vector-valued distributions. Applying 
Dencker's results, we will show that the singularities of solutions of the Dirac equation 
propagate in a simple way. This will be used to calculate the polarization set of the 
two point function for Hadamard states of the Dirac field. Our analysis makes use 



of unpublished results by Radzikowski |T8[, and our results are in agreement with 



those independently obtained by Hollands using a slightly different approach. 

This paper is organized as follows: The first chapter after the introduction gives 
an overview of the physical problems we are interested in and fixes the notation used 
throughout this work. We introduce the Klein-Gordon and Dirac fields on a curved 
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spacetime and give a definition of Hadamard states. After tliat, we come to tfie 
matliematics and review tlie metliods from microfocal analysis on which our analysis 
is based. The central definitions are the wave front and polarization sets, and the 
main results are the theorems on the propagation of singularities. As a side result, we 
will rederive the well-known properties of light propagation in a curved spacetime in 
a mathematically very elegant way. In the fourth chapter these methods are applied 
to our physical objects, and we state our main results. Finally we give an outlook 
and an overview about possible applications of our results. 

2 Quantum field theory on curved spacetime 
2.1 Notions from general relativity 

In the general theory of relativity, for an introduction see [0], spacetime is repre- 
sented by a different iable manifold {Ai,g) equipped with a Lorentzian metric g, i.e. 
a differentiable, symmetric, and non-degenerate type (0, 2) tensor field of signature 
(H ). For simplicity, we restrict ourselves to the physical four dimensions. 

As usual, the tangent bundle will be denoted TA^, and the cotangent bundle 
is T*M. For the set e T*M | ^ ^ 0} we will shortly write T*M \ 0. The 

canonical projection will always be denoted vr : T*Ai A4. Any vector bundle E 
over A4 can be pulled back to a vector bundle tt*E over T*Ai. The set of smooth 
sections of a vector bundle E over Ai will be denoted r{Ai,E), sections having 
compact support form the set Tq^MjE). 

The Levi-Civita connection is the unique connection on TAi that is symmetric 
and compatible with the metric. The associated covariant derivative is denoted V. 
The Riemann tensor is the curvature tensor of this connection which is contracted 
to the Ricci tensor and finally to the scalar curvature R. 

We will restrict ourselves to globally hyperbolic spacetimes. These are spacetimes 
that admit a Cauchy surface, that is a spacelike hypersurface H G Ai which is 
intersected by any nonextendable causal curve in Ai exactly once. In other words, 
Ai = -D(S), where we have introduced the domain of causal dependence -D(S) 
which is the set of all points x G Ai such that any causal curve through x intersects 
E. In the globally hyperbolic case the manifold is topologically of the simple form 

= M X S. Many physically interesting spacetimes like flat Minkowski spacetime, 
Schwarzschild spacetime as a black hole model, or Robertson- Walker spacetimes as 
cosmological models fall into this class. 

In order to include more general models like anti-de Sitter spacetime, this restric- 
tion can be somewhat weakened. Crucial for the following is mainly the existence 
of unique propagators as well as the existence of a spinor bundle. Since our con- 
structions are purely local, they immediately carry over to all globally hyperbolic 
submanifolds of any spacetime. 
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Globally hyperbolic spacetimes are time orientable, such that the light cone, i.e. 
the set of all nonvanishing timelike covectors in T*A4, can be separated into a forward 
and backward light cone, and V~, continuously in x. The closed light cones 
include the lightlike covectors. The time orientation induces the separation of the 
causal future J^(x) and past J~{x) of a point x e which are the sets of all points 
that can be reached from a; by a future (past) directed causal curve. 

While the tangent bundle TAi and the cotangent bundle T*Ai naturally exist for 
any spacetime A4, the existence of spinors is less trivial. In order to define spinors, 
the bundle of orthonormal frames, which is a principal fibre bundle whose structure 
group is the proper orthochronous Lorentz group, has to be lifted to a principal fibre 
bundle with the universal covering group SL(2, C) as structure group, called the 
spin structure. The Dirac spinor bundle is then defined as the associated C^-vector 
bundle. It can be shown that for any globally hyperbolic spacetime there exists 
such a spinor bundle DA4 which in this case is just a trivial bundle over Ai. The 
construction of the spinor bundle leads naturally to a covariant derivative on DAi 
and the dual bundle D*Ai which will also be denoted V. In local coordinates one 
has = dfj_ + a^, where the 4x4 matrices can be expressed in terms of the 
Christoffel symbols and the Dirac matrices. For details see ^j. 

2.2 The free Klein-Gordon field 

The simplest example of a quantum field theory is the free scalar field, satisfying the 
Klein-Gordon equation 

{Dg + m^)(j){x) = 0. (1) 

This is already the covariant field equation on a generic spacetime {Ai,g) if Og is the 
covariant wave operator Og = g^'^'W of the metric g. The positive real constant 
m plays the role of a mass. 

Because of the global hyperbolicity of our spacetime, the Cauchy problem for the 
Klein-Gordon operator has a unique solution, and there exist unique retarded and 
advanced propagators, A^et and Aq^^, such that 

{Ug + m^)Aret,adv = Aret,advi^g + "^^) = 1 (2) 

and supp{Aret,advf) *^ "^""^(supp/). We define the commutator function A = Aret — 
Aadv and identify the operator A with the distribution 

A : {f,g) ^ A{f,g) = J f{x){Ag){x)dfi{x), (3) 

where d/i is the volume element on Ai. 

In the quantum theory the field (p is represented as an operator valued distribu- 
tion such that we have smeared field operators </>(/), / G C^{Ai), acting on some 
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Hilbert space with a vacuum vector fi. They are required to satisfy the Klein-Gordon 
equation and the canonical commutation relations 

[<j>{f),<P{g)] = zA{f,g)-l. (4) 

Since we investigate a free theory, we expect that all vacuum expectation values 
of products of field operators are determined by the two point function 

A{x,y) =<n,(j){x)(j){y)n> (5) 

which is assumed to be a distribution A G T>'{Ai x A4). 

Canonical quantization of the free Klein-Gordon field can be carried through in 
close analogy to quantization on Minkowski spacetime (see for example 0), and one 
ends up with a Fock space representation of the field in terms of creation and 
annihilation operators. However, the problem is that there is no preferred mode 
expansion in terms of plain waves, because the Fourier transform is not invariant un- 
der general coordinate transformations. The crucial point is the ambiguous division 
into positive- and negative-frequency modes. Therefore this construction is highly 
non-unique, and there is a large amount of unitarily inequivalent representations, as 
was first observed by Fulling [Q]. A priori we do not know the 'real' vacuum state 
Q, and as a consequence there is no preferred interpretation of the theory in terms 
of particles. If the spacetime is asymptotically flat, one can carry over the particle 
interpretation from the flat part to the whole spacetime, but on a generic spacetime 
one has to face the fact that quantum field theory is primarily a theory of fields, and 
no unique particle interpretation can be expected. 

From an axiomatic viewpoint, the vacuum state on Minkowski spacetime is sin- 
gled out by the requirement of its Poincare invariance and the spectrum condi- 
tion which is the requirement that the energy-momentum operator (the generator 
of translations) takes its spectrum in the closed forward light cone. In the absence 
of Poincare invariance on a generic spacetime these conditions no longer make any 
sense. While one will still demand invariance under the symmetries of the spacetime 
(if there are any), there is no simple analogue to the spectrum condition since there 
are no generators of translations. Thus one central task of quantum field theory on 
curved spacetime is to characterize physical states by finding a generalized spectrum 
condition for generic spacetimes. 

Usually the construction sketched above is performed in the algebraic approach 
to quantum field theory. There the local net of observable algebras can be defined 
in a unique way, but one still has to fix a state on the algebra (which corresponds to 
the preparation of the system). Again one restricts oneself to quasifree states that 
are fixed once a two point distribution is specified, and via the GNS construction the 
state induces a Hilbert space representation of the observable algebra. For details 
on this approach we refer to 
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2.3 Hadamard states 



In order to single out a class of acceptable quantum states, one physical requirement 
is that we should be able to formulate the semi-classical Einstein equations in order 
to describe the back reaction of the quantum field on the spacetime geometry: 

Rfiu-^9fiuR = 8iT{T^^). (6) 

Here, (T^^) is the expectation value of the energy momentum tensor of the quantum 
field. Thus physical states have to allow a unique definition of this expectation value. 

It could be shown ^ that this requirement is fulfilled if we demand that the two 
point function takes the following form which was first given in and since then 
intensively studied (see [0 and the references therein): 

M^^y) = (^^^^^ + v{x,y)\n\a{x,y)\+w{x,y)^ , (7) 

where a{x, y) denotes the quadratic geodesic distance, and u, f , and w are some 
smooth functions. 

A state whose two point function is of this form is called a Hadamard state. The 
field equation and the commutation relations fix u and f , and the remaining freedom 
in the definition of a Hadamard state lies in the choice of the smooth function w. 

In order to properly define a distribution, one has to modify by an ie- 
prescription in the denominator. Also, since a is not defined globally, one must 
specify some neighbourhood in which is supposed to hold, and it is also under- 
stood that no further singularities exist for spacelike separated points. To this end, 
one demands the convergence of a series of distributions to A in a causal normal 
neighbourhood (which is basically a sufficiently small neighbourhood of a Cauchy 
surface). Thereby the Hadamard condition becomes a global condition. Although 
the choice of a Cauchy surface enters into the definition, it can be shown that the 
Hadamard condition is independent of this choice. All these details were worked out 
in p], but most of them are not essential for this work. What should be kept in 
mind is that the vacuum state on Minkowski spacetime is Hadamard, that Cauchy 
evolution respects the Hadamard singularity structure, and that two Hadamard dis- 
tributions differ only by a smooth part. 

The main lesson of this rather sketchy definition is that the Hadamard condition 
prescribes the singular short distance behaviour of the two point distribution, and it 
is obvious that it would be very useful to have at hand some powerful mathematical 
methods to characterize this singularity structure. 

2.4 The free Dirac field 

We will now briefiy outline the construction of the free Dirac field on a curved 
spacetime. For details see 0. 
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Just like the Klein-Gordon equation, we can generalize the Dirac equation to a 
spinor bundle over an arbitrary curved spacetime by replacing the derivatives and the 
Dirac matrices by their covariant counterparts. For a spinor field G T{M.^DM.) 
and a cospinor field ip G r(A^, D*M.) it reads 

{-if + m)^ = Q, = t'' V^V^; (8) 

{if + m)ij = 0, y = iy^^)^^. (9) 

Here, are the generalized Dirac matrices which satisfy the anticommutation rela- 
tions 

{7^7"} = 2^7^'^■l. (10) 

It can be shown that they form a covariant section 7 G T{M., TAi ® DAi ® D*Ai). 
Like in (|10|), we will always choose a local frame and suppress the spinor indices 
referring to D^M. ® D^Ai such that 7(x) is treated like a vector of 4 x 4 matrices. 
We can then define matrices cr^'^, 7^, 7'^7^ and 1 in the usual way, which together 
with 7^^ give a basis for DAi ® D*Ai. A further important property of 7 is that it 
is covariantly constant: V7 = 0. For any covector field ^ G r{A4,T*A4) we define ^ 
as contraction with 7, ^ = ^/j7^. 

For the quantization of the theory we again need propagators Sret and Sadv For 
globally hyperbolic spacetimes it can be shown that these propagators exist and are 
unique 0. Their difference defines the anticommutator function S = Sret — Sadv 

In the quantized theory the spinor field operator now has to be smeared with a 
smooth cospinor field of compact support in order to give an operator on the Hilbert 
space, 

ij{v), veTo{M,D*My, 

^(u), ueTo{M,DM). (11) 

The field operators are required to solve the Dirac equation and obey the canonical 
anticommutation relations, 

{2P{v)Mu)} = zS{v,u). (12) 

Again, quasi-free states are completely characterized by their two point distribu- 
tion 

u;+{x,y)=<n,ij{xmy)n> (13) 

which is now a vector-valued distribution uj^ G V{M x M,DM Kl D*M) taking 
values in the bispinor bundle DA4 Kl D*Ai. This denotes the outer tensor product: 
the fibre over [x,y) E A4 x Ai is D^Ai ® D*Ai, and the first factor in the tensor 
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product transforms like a spinor in the point x, while the second factor transforms 
like a cospinor in y. 

In the quantization procedure on Minkowski spacetime one often makes use of 
the fact that the squared Dirac operator gives back the Klein-Gordon operator. On 
a curved spacetime this is expressed by Lichnerowicz's identity: 

{-i^ + m){iyi + ni)il) = {U-^R + m'^)'ilj. (14) 

Here, R is the curvature scalar, and □ = g'^'^V^'Vi, is the spinorial wave operator 
acting on sections of DAi. Note that only the principal part of □ gives the scalar 
wave operator, and in general it contains non-diagonal terms. Nevertheless, one 
obtains the propagators Si for the Dirac equation by applying the Dirac operator to 
the propagators Aj for this spinorial wave operator, 

Si = {iy + m)Ai. (15) 

Analogously, one demands that the two point function uj^ can be extracted from an 
auxiliary two point function a), 

uj+{x,y) = {i](/^ + m)u{x,y). (16) 

Then a) is a solution of the spinorial Klein-Gordon equation, 

(^D,^y-^R + m'^u{x,y)=0, (17) 

and in order to define Hadamard states for the Dirac field one can make the same 
ansatz as in the scalar case, 

u = (- +v\n\a\+w] , (18) 



but now u, V, and w are smooth bispinor-valued functions. The rigorous definition 
was given by Kohler . Again, the singular part of the two point function is fixed 
by geometry, and the freedom in the choice of a Hadamard state lies in the choice of 
a smooth function w. 



3 Microlocal analysis 

Let us now introduce the mathematical tools we need to investigate the singularities 
of distributions. This theory of 'microlocal analysis' was developed by Hormander 
and Duistermaat in the seventies for their analysis of partial differential equations, 
see their original papers []TT], [T^ or the monographs |]TB|, . For physical applications 
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see for example [|l^, It should be noted that similar tools where independently 



developed by Bros and lagolnitzer in the context of analyticity properties of the 
S'-matrix. 

The generalization to vector-valued distributions and the definition of the polar- 
ization set was accomplished by Dencker in the eighties 
not found its way into the physical literature yet. 



1^], but it seems that is has 



3.1 Scalar distributions — the wave front set 

As test function space we will take T>{X) = C^{X), the space of infinitely differen- 
tiable functions on an open subset X C M", equipped with the usual topology. Its 
dual is the space of distributions, V'{X). 

The roughest characterization of the singularities of a distribution u G T>'{X) is 
the singular support, singsuppw. It is the set of all points x E X such that u does 
not correspond to a smooth function in any neighbourhood of x. 

But this definition can be very much refined. The following statements show that 
much information about the singularities is contained in the Fourier transform: 

1. If u is of compact support, then there exists a Fourier transform u which 
furthermore is a smooth function. 

2. A distribution u of compact support corresponds to a smooth function if and 
only if u decays faster than any power, that is if for every m G N there is a 
constant Cm G M such that 

<Cm(i + iei)"'" veGM^ (19) 



Thus in order to investigate the singularities of a distribution u G V{X) locally 
at the point x G X, we are led to analyze the Fourier transform (pu for test functions 
G C^{X) with 7^ which cut off the singularities far away from x. Then 
X G sing supp u iff there is no such that (pu falls off rapidly in all directions. 

This immediately suggests two refinements of the singular support: we may ask 

• If (pu does not decay rapidly in all directions, which are those directions in 
Fourier space which are responsible for the singularity? 

• If (f)u does not fall off rapidly, how fast does it decay? 

The first question leads directly to the notion of the wave front set, while the second 
question motivates the definition of local Sobolev spaces H^{x). In combination one 
can define if'^-wave front sets, but in the following we will only be interested in the 
C°°-case: 
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Definition 1. Let u e V'{X), X C M". A point {x,^) e X x {W \ {0}) is called 
a regular directed point of u if there is a function (f) e C^{X), not vanishing at x, 
such that for any m G N there is a constant Cm G K with 

\Ma\<Cmii + \e\r"' (20) 

for all ^' in a conic neighbourhood F C \ {0} of ^. (A neighbourhood F is called 
conic if with ^ G F all points t ■ C,, t & M+, are contained in F.) 

The wave front set of u, denoted WF{u), is the complement in X x (M" \ {0}) of 
the set of regular directed points of u. 

Example. The wave front set of the 5-distribution 6 G 'D'(M") is 

W(5) = {(0,OkGK"\{0}}. (21) 
This can directly be seen from its Fourier transform 

m) = 5(e-<«'>0) = 0(0) (22) 

which does not decay rapidly in any direction. 

Our second example shows that the wave front set does not always contain the 
whole Fourier space for all points of the singular support: 

Example. The distribution u G V(M.), u{x) = has got the wave front set 

W(w) = {(0,OKeR+\{0}} (23) 
which can be seen from its Fourier transform 

u{0 = -^v^^(0- (24) 
We will now collect some basic properties of the wave front set [|ll[]: 



Theorem 1. 1. WF{u) is a closed conic subset of X x (R" \ {0}). 

2. The wave front set is a refinement of the singular support: 

n{WF{u)) = singsuppw. (25) 
In particular, the wave front set of a smooth function is empty. 

3. For any partial differential ( or more generally any pseudo-differential) operator 
P we have the pseudolocal property 

WF{Pu) C WF{u). (26) 
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4- Let U,V G M", u G T>'{V), and x '■ U ^ V a diffeomorphism such that 
X*u G V'{U) is the distribution pulled back by x- X*u{f) = u{f o x"^). Then 

WF{x*u) = x*WF{u) = {{x-\x), yix-\xm) I e WF{u)} . (27) 

In particular, under coordinate transformations the elements of the wave front 
set transform like covectors. Thus the wave front set can be defined for distri- 
butions on differentiable manifolds by gluing together wave front sets over the 
coordinate patches, and for u G V\M.) we have WF{u) C T*M. \ 0. 

The wave front set can be defined in a second way which will be important 
for the investigation of solutions of partial differential equations as well as for the 
generalization to vector-valued distributions. But first we need some notions from 
the theory of partial differential operators (PDOs): 

A PDO A is a polynomial of partial derivatives, A = P{x, d). In momentum space 
it acts like multiplication by the same polynomial of momenta P{x,iC,) = (Ta{x,C,) 
which is called the symbol of A. The leading order a(x, ^) of a a is called the principal 
symbol ap{A). The generalization to a manifold A4 via local charts is obvious, and 
it should be noted that the principal symbol is then a well-defined function on T*Ai. 

Usually this is generalized to pseudo-differential operators (\E'DOs), where 
o"a(x,^) may be a more general but still well-behaved function, not necessarily a 



polynomial (for details see |]T4|)- In this work we will only deal with classical \E'DOs 
whose symbol is an asymptotic sum of homogeneous terms. The set of classical 
\E'DOs on ^A of order m will be denoted L^{A4). 

The characteristic set of a \E'DO A G L"^{Ai) with principal symbol a = ap{A) 
is defined as 

char(A) = { (s, G T*M \ | a{x, = O} (28) 

and can be interpreted as the set of all directions ^ suppressed by A to leading order, 
at a point x. 

Now the following theorem makes precise the intuitive statement that the sin- 
gular directions of a distribution (which make up the wave front set) are just those 
directions that have to be suppressed to leading order by any operator that maps 
the distribution to a smooth function. 

Theorem 2. For the wave front set of a distribution u G V'{M.) we have 

WF{u) = Pi char(v4), (29) 

A«gC°=(A4) 

where the intersection is taken over all pseudo-differential operators A G L°(A^) with 
Au G C°°{M). 



11 



This property of the wave front set is particularly useful if we know such an 
operator A with Au G C°°, especially if u is a solution of a partial differential equation. 
The first part of the following theorem is then obvious, but for certain operators one 
can make an even stronger statement which goes under the name 'propagation of 
singularities theorem' and was first proved by Duistermaat and Hormander [Q: 

Theorem 3. Let P G L™(A^) be a ^DO on M. with principal symbol p. If u & 
V'{M) such that PueC°°{M), then 

WF{u) C char(P). (30) 

// furthermore P is of real principal type, then WF{u) is invariant under the flow 
generated by the Hamiltonian vector field of p. 

Here we have imposed a natural restriction on P in order to obtain a real Hamil- 
tonian and a non-degenerate Hamiltonian flow: 

Definition 2. A pseudo-differential operator P G L"^{Ai) is said to be of real prin- 
cipal type if its principal symbol p{x, ^ ) is real and for p = the Hamiltonian vector 
field Hp, 

Hp{x, = Y1 O-^xM - 9^mp(x, 0%, (31) 

does not vanish nor does it have the radial direction, that is Hp ^ 

Thus the wave front set of a distribution u with Pu = is made up of integral 
curves of Hp in char(P) which are also called the null bicharacteristics of P. Their 
projections onto M. are called the bicharacteristic curves of P and constitute the 
singular support of u. 

To illustrate this theorem, let us take a look at a physical example, the wave 
operator on a curved spacetime: 

Example. Let {M,g) be a spacetime, /, 6 G C°°(A<), such that ^ f{x) G M Vx G 
M, and a G r{M,TM). The operator 

P = /□ + a^^V, + b, □ = g^-'V^V,, (32) 

has got the principal symbol p{x,C,) = — f {x)g'^^ {x)S,^iu = —f ■ i"^- Hence, P is of 
real principal type. 

For (x, ^) G char(P) we see that the covector ^ has to be lightlike. Furthermore, 
for = the function p(x, ^ ) = — is a well-known Hamiltonian that generates 
the null geodesies which is not affected by the nonvanishing factor /. Therefore the 
bicharacteristic curves of P are null geodesies x(r). 
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From Hamilton's equations we obtain 

dr d^^ 



-2f{x{r))g^^i. (33) 



which means that ^ is tangent to the geodesic. By reparameterization of the geodesic 
one can achieve ^(r) = ^(t) for all values of the curve parameter r. The factor, 
however, is irrelevant in what follows since the wave front set is conic in ^. 

Thus the null bicharacteristics of P are curves (x, ^)(r) G T*M. \0 such that x(r) 
describes a lightlike geodesic and ^(r) is tangent to x(r) for all values of r. 

Applying theorem |^, we conclude that the singular support of a distribution u 
which solves the wave equation, i.e. Pu = 0, is a union of null geodesies, while the 
wave front set additionally indicates the covectors tangent to the geodesies in the 
singular support. Note that the geodesies are directed and that the wave front set 
may contain only the future or past directed covectors. 

If we take light rays as distributional solutions of the wave equation, this result can 
be interpreted as the well-known fact that light propagates along lightlike geodesies. 

3.2 Vector-valued distributions — the polarization set 

We will now go over to vector- valued distributions u G T''(A^,C^), i.e. vectors 
u = {ui)i=i...N of distributions Ui G V'{A4). This section is mainly a review of . 



The wave front set of a vector- valued distribution u = (m,) G V'{Ml, C^) is just 
defined as the union of the wave front sets of all its components: 

WF{u) = [jWF{ui). (34) 

i=l 

Because of theorem ^ this is nothing but 

N 

WF{u) = [j Pi char(A). (35) 

j=l Aui€C°°{M) 

The wave front set does not contain any information about the components of 
the distribution that are singular. In order to specify the singular directions in the 
vector space C^, one could consider vector- valued operators that map the vector- 
valued distribution to a smooth scalar function, instead of just looking at scalar 
operators mapping the individual components to smooth functions. This approach 
leads to the definition of the polarization set: 

Definition 3. The polarization set of a distribution u G T>'{A4, C^) is defined as 

WFpoi{u)= fl Ma, (36) 

Au£C°°{M) 
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Na = {{x,^;w) e {T*M\0) X C^|w e kera(x,0}, (37) 

where the intersection is taken over all 1 x systems A G L^{Ai)^ of classical ^DOs 
with principal symbol a, and ker a(x, ^) is the kernel of the matrix a{x, 

Obviously, {T*M \ 0) x {0} C WFpoi{u) for any u e V'{M,C^). Furthermore, 
the polarization set is closed, linear in the fibre and conic in the ^-variable. 

For scalar distributions (A^ = 1) the polarization set contains the same infor- 
mation as the wave front set. For arbitrary N, we get back the wave front set by 
projecting the nontrivial points onto the cotangent bundle: 

Theorem 4. Let u e V'{M,C^) and 7ri,2 : T*M x ^ T*M be the projection 
onto the cotangent bundle: tti 2(2;, w) = (x,^). Then 

TTiAWFpoiiu) \ {T*M X {0})) = WF{u). (38) 

In this way the polarization set is a refinement of the wave front set. In addition, 
it contains information about the directions in the additional vector space in which 
the distribution is singular. 

Example. Let u = (^1,1*2) G V^M^C^) and {y,rj) ^ WF{ui). Then 

WFp,i{u) C {(x, ^; (0, z)) e (T*M \ 0) x C^ z e C} (39) 

over a conic neighbourhood of {y,i]). 

The polarization set indicates only the most singular directions, even if the projec- 
tion of the distribution on other directions is also singular, as the following example 
shows: 

Example. Let u = {v,Av) G V'{W,C^), where v e V'{W) and A is the Laplacian 
on W. Then 

WF,a{u) C {{x,i;{0,z)) e (r*M" \ 0) x C^, ^ e C}, (40) 

because Aui — U2 — O and ker(7p(A, — l)(x, ^) = ker(— ^^,0) = {(0,2;), z e C}. 

Our third example shows that the direction of the strongest singularities in C''^ 
can depend on the singular direction in the cotangent bundle: 

Example. Let u = V(5(2) = {d^,5^'^\ d^J^'^^) e ©'(M^C^). Then 

WF,,i{u) c { (0, e; A • e {T*R^ \ 0) X C^ A e C}, (41) 
since u solves the equation (—5x2, ■ u = 0, and ker(— i^2, ^Ci) = ^ " (^1)^2) for 

The polarization set has got an interesting transformation property under the 
action of a system of partial differential operators: 
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Theorem 5. Let A be an M x N system of pseudo-differential operators on M. with 
principal symbol a{x,^), and u G P'(A^,C^). Then 

a{WFp,i{u)) C WFj,,i{Au), (42) 

where a acts on the fibre: a{x,^;w) = {x,C,;a{x,C,)w). 

If E is an N X N system of pseudo-differential operators on A4 and its principal 
symbol e{x,C,) is not characteristic at {y,r]) G T*A4 \ 0, i.e. det e{y,ri) ^ 0, then 

e{WFpM) = WFp,i{Eu) (43) 

over a conic neighbourhood of {y, rj). 

Note the different order of the inclusion in (|4^ ) compared to (^61). Since the polar- 
ization set only indicates the most singular directions, the projection on orthogonal 
directions can change the polarization set substantially. 

From (^3]) we learn that the polarization set behaves covariantly under a change 
of basis in the vector space C^. Thus the polarization set of a distribution u G 
T>'{/iA,E) taking values in a vector bundle E ^ M can be defined by gluing to- 
gether the polarization sets of local trivializations. This gives a well-defined subset 
WFpoi{u) C 71* E of the vector bundle E lifted over the cotangent bundle. 

Like in the scalar case, we are interested in the propagation of singularities for a 
solution of a system of partial differential equations. Again it turns out that there 
is a powerful theorem for a restricted class of operators. 

Definition 4. An N x N system P of pseudo-differential operators on A4 with 
principal symbol p{x,C,) is said to be of real principal type at {y,ri) G T*Ai \ if 
there is an x symbol p{x,C,) such that 

P{x,^)p{x,^)=qix,^)-1 (44) 

in a neighbourhood of {y,ri), with a scalar symbol q{x,^) of real principal type. We 
say that P is of real principal type in (9 C T*A4 \ if it is at all points {y, rj) G O. 

The choice of p and q is not unique, but it can be shown that the only freedom 
is the multiplication by a smooth non-vanishing function on T*Ai \ 0. 

From now on let P be an A^ x A^ system of classical pseudo-differential operators 
P of real principal type of order m, u a. solution of Pu = 0, and locally p{x, ^) and 
q{x,C,) are chosen as above such that 

pp = q ■ 1. (45) 

We define the set 

= [{x,0 e T*A^ \ 0| detp(x,0 = O}. (46) 
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Because of (|45|) , we have Qp = q^^{0) locally. If Pu G C°°, for the existence of 
nontrivial elements {x,^]w) G A/p in the polarization set over a point (x,^), by 
definition it is necessary that w G kerp(x, that is, kerp(x, ^) 7^ {0}. This means 
that (x,^) and thus the whole set WF{u) has to be a subset of Qp. 

It can now easily be shown that the wave front set is a union of null bicharac- 
teristics of q in Qp. Though the bicharacteristics are not independent of the choice 
of q, we will just call them the bicharacteristics in without referring to a special 
choice of q. A different choice of q only leads to a multiplication of the ^-variable by 
a scalar factor which is irrelevant since the wave front set is conic in ^. In fact, in 
general it is not even possible to choose q globally, and several coordinate patches 
must be glued together anyhow. 

Once we know the wave front set of the distribution u, the remaining task is to 
calculate the polarization vectors over the points in the wave front set. It turns out 
that these vectors follow a simple parallel transport law along the bicharacteristics 
that form the wave front set. 

This parallel transport, which we will introduce in a second, does not only depend 
on the principal symbol as in the scalar case, but also the subprincipal symbol 
contributes. The symbol of P is a sum of homogeneous terms, 

(^iP)ix,0 =p(x,0 +Pm-i(a;,0 +Pm~2(a;,0 + (47) 

where p = crp{P) is the principal symbol, and pj is homogeneous of order j. The 
subprincipal symbol is now defined as 

Furthermore let 

_ dp{x, dp{x, ^ dp{x, dp{x, Q 



^ d^f^ dxf^ ^ dx^ d^^ 



be the Poisson bracket and 



the Hamiltonian vector field of q. 

Definition 5. For a smooth section w of the vector bundle {T*M. \ 0) x C^, we 
define Dencker's connection as 

Dpw = HqW + ]:{p,p}w + ipp^w. (51) 
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This is a partial connection along the bicharacteristics in fip, that is a connection 
on all C^-vector bundles over these bicharacteristics. Furthermore, it can be shown 
that (51) defines a partial connection in A/p, this means that for every vector field w 



along a bicharacteristic 7 in fip we have Dpw G kerp along 7 if and only if w G kerp 
along 7. The equation Dpw = can then be solved with (x, ^; w) G A/p, which is a 
necessary condition for elements of the polarization set. 

Again, the definition of the connection depends on the choice of p and q, but it 
can be shown that a different choice changes a solution of the equation Dpw = in 
A/p only by a scalar factor. Therefore, we define a Hamilton orbit of a system P of 
real principal type as a line bundle L C A/p|-y over a bicharacteristic 7 in f2p which 
is spanned by a section w satisfying the equation Dpw = 0, i.e., that is parallel with 
respect to Dencker's connection. These Hamilton orbits are then independent of the 
choice of p. 



Now we are prepared to state the main result of |jT^, the theorem on the propa- 
gation of singularities for vector-valued distributions: 

Theorem 6. Let P be an N x N system of classical pseudo-differential operators 
over a manifold M., and u G 'D'{M.,<C^). Furthermore, let P he of real principal 
type at {y,r]) G flp, and {y,r]) ^ WF{Pu). Then, over a neighbourhood of {y,ri) in 
flp, WFpoi{u) is a union of Hamilton orbits of P. 

Though it is not obvious from the definition (|5T]), in interesting cases Dencker's 
connection takes on a very simple form. We will illustrate this by a physical example 
that was already discussed by Radzikowski [|18| in a slightly modified way. 

Example. We investigate Maxwell's equations on a spacetime {/iA,g). For a vector 
field A G r(A4, T/A) in Lorentz gauge (Vp^A'^ = 0), in a local coordinate frame they 
read 

ngA-" - R\A^ = 0, = g'^'VpV,. (52) 

We want to calculate Dencker's connection for the operator 

P\ = n/^ - R\ = g'''{6;dpd^ + T%pd„ - 5;r\.9A + T\^dp) + . . . , (53) 

where the dots stand for lower order terms that do not contribute to Dencker's 
connection. 

The principal symbol of P is 

p{x.i)\ = -g"'{x)UJ';. (54) 

so that we can choose 

P\ = v^5', q{x, = -V^g'^U. (55) 
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to get pp = q-1. The factor is introduced for no obvious reason, but it turns out 
that it simphfies the calculation. We see that g is a scalar symbol of real principal 
type, therefore P is of real principal type. 

The bicharacteristic curves of q are again the null geodesies, and the space Qp = 
g~^(0) consists of all lightlike vectors in T*A4 \ 0. Since p oc 1, the space kerp(a;, ^) 
on which Dencker's connection acts is the whole tangent space. 



Putting together the different terms in (|5lD and using well-known identities for 



the Christoffel symbols and Hamilton's equations, we arrive at 
{DpwY = ^HqW + ^{p,p}w + ipp^'iu^ 

= ((rVr)vJr, (56) 

for a vector field w{x{t),C,{t)) along the null bicharacteristics. That is, Dencker's 
connection for P (with respect to our choice of p) is just the usual Levi-Civita con- 
nection, lifted over the cotangent bundle and restricted to the bicharacteristics. Thus 
a vector field over one of the bicharacteristics is parallel with respect to Dencker's 
connection (and thereby generates a Hamilton orbit) iff the projected vector field 
over the characteristic curve is parallel with respect to the Levi-Civita connection. 

Thus, following the theorem on the propagation of singularities, the polarization 
set of a solution of Maxwell's equations is a union of such Hamilton orbits, that is it 
consists of curves (x,^; w){t) G 7r*TA^ such that 

• x{t) describes a null geodesic, 

• ^(r) is tangent to the geodesic and 

• w{t) is parallel transported along the geodesic. 

In physical language, we have reproduced the well-known result on the propaga- 
tion of light in a curved spacetime: light travels along lightlike geodesies, while the 
polarization vector is parallel transported along the path. 

In particular this means that spacetime curvature cannot lead to double refraction 
which is only possible if matter effects are introduced into Maxwell's equations. 
Double refraction may then occur at points where the operator P is no longer of real 



principal type. This was also investigated by Dencker [IS 



4 Singularity structure of the two point function 

We are now going to investigate the singularity structure of the two point function 
of Hadamard states using the microlocal techniques that were introduced in the 
previous section. 
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The wave front set of a scalar Hadamard distribution was first calculated by 



Radzikowski ||2^. After we have reviewed his results, we will derive an analogous 
result for the wave front set of the two point function of the free Dirac field and 
furthermore calculate its polarization set. 

4.1 Klein-Gordon field 

Let us first state the main result by Radzikowski [2 



Theorem 7. A quasi-free state of the Klein- Gordon field on a globally hyperbolic 
spacetime M. is a Hadamard state if and only if its two point distribution A has got 
the following wave front set: 

WF{K) = {ix,y;^, -v) E T\M x A^) \ | (x, ~ (z/, r/), ^ e F^}. (57) 

The equivalence relation (x,^) ~ {y,v) nT'^ans that there is a lightlike geodesic 7 
connecting x and y, such that at the point x the covector ^ is tangent to 7 and rj is 
the vector parallel transported along the curve •y at y which is again tangent to 7. 
On the diagonal {x,C,) ~ {x,ri) if ^ is lightlike and C, = rj. 

Thus the two point distribution of a Hadamard state is singular only for pairs 
of points that can be connected by a lightlike geodesic. In addition, the condition 
^ G expresses the fact that only positive frequencies contribute and can be viewed 
as the microfocal remnant of the spectrum condition. The obvious advantage of our 
formalism from microlocal analysis is that the spectrum condition (for the free scalar 
field) is now formulated also for curved spacetimes. 

This microlocal characterization of the Hadamard condition opens the door for 
further investigations of Hadamard states using the powerful tools from microlocal 
analysis. As a first example, in |2l| it was shown that certain states that were known 



long before, the ground- and KMS-states on static spacetimes and adiabatic states 
of infinite order on Robertson- Walker spacetimes, are indeed Hadamard states. 
We will now briefly review Kohler's proof |T0| that the two point function of 



a Hadamard state has got the given wave front set, since it will lead us in the 
generalization to the Dirac field. 

First we take a look at the two point function on fiat spacetime: 

Theorem 8. On Minkowski spacetime the two point function of the free Klein- 
Gordon field in the vacuum state has got the wave front set 

WF{A) = {{x,y-^,-OeT*{R'''xR''^)\xj^y,{x-yY = 0, 

eii(x-z/),eo>o} 

U {{x, X- e, -0 G T*(Ri'3 X Ri'3) I ^2 ^ 0, eo > 0}, (58) 
i.e. it is of the form (j3^. 
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A complete proof can be found in |jT5|, or the wave front set can be read off from 
the well-known Fourier transform 



A(e, v) = (27r)-^5(e + vWi^one - m'). (59) 

The extension of this theorem to arbitrary spacetimes is done by deforming the 
spacetime such that the metric is flat in one part of the spacetime. The known sin- 
gularity structure in the flat part can then be shifted to the curved part by applying 
Hormander's theorem on the propagation of singularities. 



The proof makes use of the following theorem which was proved in [|lOl based on 
ideas from ||22|| : 

Theorem 9. Let {Ai,g) be a globally hyperbolic spacetime and x a point on a Cauchy 
surface (Z M.. Then there is a neighbourhood U of x and a globally hyperbolic 
spacetime {Ai,g) with the following properties: 

1. A causal normal neighbourhood M of H with U d M is isometric to a neigh- 
bourhood M in M., and the isometry p : J\f ^ J\f maps S to a Cauchy surface 
tcAf. 

2. There is a spacelike hypersurface S together with a neighbourhood U in M. such 
that g, restricted to U, is the Minkowski metric and p{U) = U G -D(S). 

The two point function A which is given on the neighbourhood Af of the Cauchy 
surface S and thereby also on = p{Af) now induces a Hadamard two point distri- 
bution A on the whole deformed spacetime M. such that A\_\fxAf = P*{Mmxm)- 

In the flat part of A4 we already know the wave front set of the distribution A. 
Since only local properties of the spacetime in a neighbourhood of a Cauchy surface 
enter into the Hadamard condition and the metric on U is flat, on f/ x f/ we have 
because of theorem 1^: 

ix,y;i,-f,)eWFiA\^^^) ^ (x, i) iy,v), ^ eVt . (60) 

A solves the Klein-Gordon equation in the second variable, therefore we can apply 
theorem ^ on the propagation of singularities for the operator 1 ® (□ -t- m^) to get 
the singular directions of A for points in U x U. From our investigation of the wave 
operator in the previous section, we learned that two points (x,^) and (y, r/) lie on 
the same null bicharacteristic for the operator □ + iff {x,C,) ~ {y,v)- Because of 
U C D{U), every null geodesic through U intersects U, and theorem |^ tells us that 

{x,y;i-r])eWF{Ay^^) ^ {x,y; i -fj) e WF{A\(,^^), {y,r]) ^ {y,fi) 

^ {y,v)-ix,i),ievt (61) 
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The same argument for the operator (□ + m^) (g) 1 gives us the wave front set for 
points in U X U: 

{x,y;^,-7^)eWF{A\^^^) ^ (x, ~ (x, 6, (z/, ^) ~ 0, ^ e Ft, (62) 
and therefore 

WF{A\^,^) = {ix,y;^, -r/) G T*{U x f/) \ | ~ {y,v),^e F^}. (63) 

Now A is the distribution pulled back by the isometry p, and because of theo- 
rem [^.4 also the wave front set of A in f/ x [/ is of this form. 

Using the same line of arguments again in the undeformed spacetime Ai, we 
finally obtain the wave front set of A for arbitrary pairs of points {x,y) & A4 x Ai: 

WF{A) = {{x,y;^,-r])e T\M x A< ) \ | (a;, ~ (y, ^) , ^ e F+ } , (64) 
and we have finished the proof. □ 

4.2 Dirac field 

Let us now turn to Hadamard states of the free Dirac field. We first state our main 
result: 

Theorem 10. Let u be a Hadamard state of the free Dirac field on a globally hyper- 
bolic spacetime M. . Then its two point function uo^ has got the following wave front 
and polarization sets: 

WF{u;+) = {{x,y-^,-r])eT\MxM)\0\{x,O-^{y,v),^eVt}, (65) 

WFpoiiuj-^) = {ix,y;^,7];w)en*{DMMD*M)\ 

(x, y; e, v) e WF{uj+y, (1 ® J,{x, y))w = A ■ jf, A G C}. (66) 

Here, J^[x^y^ : D*Ai denotes the parallel transport in D*A4 along the 

geodesic 7 connecting x and y, such that ^ is tangent to 7 in the point x. 

The proof is similar to that for the scalar case: We will first calculate the polar- 
ization set on Minkowski spacetime and afterwards shift the result to our spacetime 
Ai by making use of Dencker's theorem on the propagation of singularities. 

Theorem 11. The two point function uj~^ of a Hadamard state of the free Dirac field 
on Minkowski spacetime Ai = M^'^ is (up to a smooth part) of the form 

a;+ = [(i^ + m)®l](A-l), (67) 

where A is the two point function of a Hadamard state of the free Klein-Gordon field. 
Rs polarization set is 

WFp,i{u+) = 

{ (x, 2/; e, r^; A • ^ G 7T*{DM M D*M) \ (x, y; e, v) e WF(A), A G C}. (68) 
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Proof. The bundle DAi Kl D*Ai over Minkowski spacetime Ai = M^'^ is the trivial 
bundle [Ai x At) x C^^'*. Thus we can simplify our notation by identifying the 
spinor spaces over all points. Also, the covariant derivative coincides with the partial 
derivative of the individual components, and the curvature scalar vanishes. Therefore 
the functions u and w in the definition (0) of a Hadamard distribution have the 
form u = u-l, v = v- l with the corresponding scalar functions u and v. That 
is, the auxiliary two point function of any Hadamard state of the free Dirac field on 
Minkowski spacetime is (up to a smooth part) a multiple of the unit matrix, and the 
nonvanishing components are two point functions of scalar Hadamard states. Since 
these are fixed up to a smooth part, the auxiliary two point function is of the form 
a) = A ■ 1 with a scalar Hadamard distribution A (up to a smooth part), and we have 

= [(i^ + m) ® 1](A- 1). (69) 

From theorem |^ we know that the wave front set of A ■ 1 has the form (pTf), and 
its polarization set is obviously 

WFpoiiA ■ 1) = { (x, y; r^; A ■ 1) | (x, y; ^, r]) e WF{A), A G C}. (70) 

We obtain from A ■ 1 by application of the operator 

A=[{i^ + m)^l] (71) 

with principal symbol a{x, y;C,,i]) = —^®l. Following theorem |1|, this does not en- 
large the wave front set. In addition, from theorem ^ we have the following restriction 
on the polarization set of uj~^ = y4(A ■ 1): 

WFpoi{uj^) ^ {{x,y;^,r];a{x,y;^,r])w)\{x,y;^,r];w) eWFpoi{A-l)} 

= {{x, y- {x, y- e WF{A), A G C}. (72) 

Now ^ 7^ if {x, y; ^, rj) G WF{A), and since the projection of the nontrivial part of 
the polarization set onto the cotangent bundle gives back the wave front set, we see 
that by the action of A the wave front set does not become smaller. Thus uj~^ has 
got the same wave front set as the scalar Hadamard distribution A. 

The equality for the polarization set in (0), however, does not follow from the- 
orem 1^, since the operator A is characteristic just in the interesting points: we have 
deta(x,2/;e,r/) = -{e? = for {x,y;^,r]) G WF{A). 

Instead, we can give a direct calculation of the polarization vectors w over a point 
{x,y]^,ri) in the wave front set. Because of (|69|), such vectors have to be a linear 
combination of the unit and gamma matrices, and we can set w = al + jS^'j". 

By definition, if a is the principal symbol of an operator A with Auj~^ = 0, a point 
(x, y; C,, 1]] w) can only be contained in the polarization set of uj~^ if a{x, y; ^, ri)w = 0. 
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We know such an operator, since the two point distribution is a solution of the Dirac 
equation 

(-^^, + mV+(x,y) = 0, (73) 

and we obtain 

^=^-w = i,Y{al + M") = ai,Y + i,pn - I Y^i^^Pu - ^uPf.)(T^'', (74) 

where a^'^ = ■^['y^,l'^]- Hence, since ^ 7^ 0, we must have a = and /3 = A ■ ^ with 
A e C. This means that the polarization vectors over points {x,y]C,,Ti) G WF{u)~^) 
are indeed only the vectors of the form w = X ■ ^. 

This completes the proof. □ 

We are now going to give a generalization to arbitrary globally hyperbolic space- 
times by using Dencker's theorem on the propagation of singularities of vector-valued 
distributions. First we calculate Dencker's connection for the Dirac equation: 

Theorem 12. Dencker's connection Dp for the Dirac operator P = {—i'^ + m) 
acting on the Dirac bundle DM. over an arbitrary spacetime M. can be chosen to 
coincide with the spin connection lifted to n*DAi and restricted to the null geodesies. 
The analogous statement holds for the Dirac operator on the dual bundle. In partic- 
ular, the connection is independent of the mass m. 

Proof. A first (rather involved) proof for the Dirac equation for spinor half densities 
was given by Radzikowski |jl8| . Here we give a simple proof that holds for the Dirac 
equation on DA4. 

Once we haven chosen coordinate frames for the tangent and cotangent bundle, 
the operator 

+ m = -iY{d^, + af,)+ml (75) 
has principal and subprincipal symbols 

P{x, = ^^.l^ix), p'{x, = -i7^(x)or^(x) + ml - ^^d^l^ix). (76) 
We choose 

Pix, = \/^gi^^,rix), qix, = \/^gix)g''''ix)^,^,, (77) 

and see that pp = q ■ 1, thereby the operator is of real principal type. The charac- 
teristic curves of q are again the null geodesies. Dencker's connection acts on spinor 
fields w along the bicharacteristics that lie in the kernel of p{x,C,) for every point 
(x,^), that is they satisfy ^w{x,^) = 0. 



23 



Putting together the different terms in Dencker's connection and using = 0, 
Hamilton's equations and the properties of the Christoffel symbols and the Dirac 
matrices, yields 

Dpw = HgW + ^{p,p} + ipp" 

^ + x^cr^^ w = in*Vr)w, (78) 

which is just the covariant derivative of w (as a vector field along the geodesic) along 
the geodesic line with respect to the spin connection. 

The proof for the adjoint operator is done in the same way. □ 

The proof of theorem [1^ now proceeds like in the scalar case. The spacetime is 
deformed such that one part U carries the Minkowski metric. By the isometry p, the 
undeformed part of the deformed spacetime inherits a spin structure which can be 
extended to the whole spacetime, since in the globally hyperbolic case all spin struc- 
tures are trivial bundles. Because the metric stays unchanged in a neighbourhood of 
a Cauchy surface, the two point function of the Hadamard state on the undeformed 
spacetime induces a Hadamard two point function on the whole deformed spacetime. 

The polarization set in the flat part U was calculated above. On the diagonal we 
have 

(x,x;e,r);«;)eW^FpoK^+lc/x[/) ^ -fj = ^ EVt = 0;w = X ■ I X e C. (79) 

In order to calculate the polarization vectors w over a point {x,x) E U xU in the 
curved part of the spacetime, we make use of the fact that the two point distribution 
solves the equation 

Pcj+ = {{-iy + m) (g) 1 + 1 (g) (iy + m))u;+ = 0. (80) 

The curves (x, x; ^, — ^)(r) such that x(r) is a null geodesic and ^(r) is tangent to 
the curve x{t) for any r are null bicharacteristics of this operator. As a corollary to 
theorem |12|, one can easily see that Dencker's connection for P is Dp = 7r*(VT ® Vr) 
along these curves. 

According to theorem the polarization set of u!~^ consists of Hamilton orbits for 
the operator P. These are sections w of the bundle TT*{DAi Kl D*Ai) along the null 
bicharacteristics that are parallel with respect to the connection Dp. In our case, 
for the bicharacteristics {x, x, ^, — ^) this just means that the polarization vectors 
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w = X ■ ^ are parallel transported along the geodesic curves {x, x){t). Now we have 
V ■ ^ = along the geodesies, since ^ is covariantly constant (the tangent vector of a 
geodesic is parallel transported along the curve) as well as the Dirac matrices. Thus 
the polarization vectors retain their form w = X - ^, and the polarization set over the 
diagonal in the curved part is of the same form as in the flat part of the spacetime 
and consists of points (x, x;^, — ^; A ■ ^) such that ^ is lightlike and in the forward 
light cone, with arbitrary A G C. 

Because of theorem ^ this form of the polarization set is conserved when the two 
point distribution is pulled back to the undeformed spacetime via the isometry p. 

Finally, we obtain the polarization vectors away from the diagonal by shifting 
the second spinor to the second point: for the operator P = 1 ® (i^ + m) Dencker's 
connection is simply Dp = 1 ® V. Thus a point (x, y; ^, — r/; w) is contained in the 
polarization set of uj^ if and only if the parallel transport along the null geodesic 
connecting x and y, such that ^ and rj are tangent to the geodesic, shifts a polarization 
vector A • ^ over the point (x, x, ^, — ^) to w. 

This is just the statement of theorem 0. □ 

To summarize, we have proved that, like the wave front set, the polarization set 
of the two point function of a Hadamard state is uniquely fixed by the underlying 
geometry, and that the fibre over each point of the wave front set is one-dimensional. 
Thus the polarization set takes on the minimal possible form. 

It seems now natural to define a Hadamard state by the polarization set of its 
two point distribution. A first step in this direction was undertaken by Hollands |]23 



who takes equation (q^) for the wave front set, which is a result in our approach, as 
the definition of a Hadamard state for the free Dirac field. It is then shown that this 
condition on the wave front set already fixes the polarization set. 



5 Outlook and conclusion 

In this paper we have demonstrated the usefulness of a microlocal characterization 
of the distributions that are of interest in quantum field theory on curved spacetime. 
However, this is not confined to the investigation of quasifree Hadamard states. 

In the naive perturbative construction of interacting quantum field theories one 
encounters formal products of distributions that are a priori not well-defined. Over 
the years, one has learned how to deal with the divergencies in these formal expres- 
sions by different renormalization techniques, but the understanding of the singular- 
ity structure of the distributions which are involved can lead to some further insight 
where these divergencies originate. 

For example, in scalar 0^-theory the perturbative construction of the S'-matrix 
involves the time-ordered two point distribution, whose formal expansion in terms 
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of Wick products of free fields contains terms like 



T2{x,y) = ci{i/\F{x,y)f ■.<P{x)<P{y): +C2{i/^F{x,y)f ■.(t)\x)(t)\y): + . . . , (81) 
where 

AF = iK + A+ (82) 

is the Feynman propagator. 

Because of the singularities of distributions, it is impossible to define a reasonable 
product on the whole space of distributions, especially the products of Feynman 
propagators in (^) do not exist as well-defined distributions. This is the reason for 
the well-known ultraviolet divergencies in these expressions. 

For certain well-behaved distributions, however, it is possible to define a product, 
and it can be shown that a simple condition on the wave front sets of the factors is 



sufficient for its existence |Tl|: If the singular directions in the wave front sets of the 
two factors over the same point do not add up to zero, the product can be defined. 

The wave front set of the Feynman propagator Aj? of a Hadamard state has also 
been calculated by Radzikowski [PD]], 

WF{^f) = {(x,y;e,-r7)GT*(A^xA^)\0| (x,eH(y,r/),a;^y, 

U {(x, X- e, -0 G T\M X M) \ 0}. (83) 

The condition ^ G if x G J^{y) now ensures the existence of the products in (^Tp 
at all points away from the diagonal, while over the points of the diagonal the wave 
front sets of the factors are too large. Therefore the only ambiguity in the definition 
of the product lies in the extension of the resulting distribution to the diagonal, 
and renormalization can be viewed as the process of extending the time-ordered 
distributions to the whole spacetime. 

This approach to renormalization theory has not only the advantage of being 
mathematically elegant and rigorous, but it also works entirely in configuration space. 
Thus it can be extended to a generic spacetime, where all the powerful renormal- 
ization techniques in momentum space cannot be applied because of the lack of a 



global momentum space. Indeed, in [|2J] it could be shown that scalar field theories on 
curved spacetimes can be renormalized under the same conditions as on Minkowski 
spacetime. 

As the main new result of this work, we have calculated the polarization set for 
the two point function of a Hadamard state of the free Dirac field. The Feynman 
propagator can be investigated in the same way, and on Minkowski spacetime it can 
easily be shown that the polarization vectors are of the same form as for the two point 
function. The generalization to a curved spacetime, however, is complicated by the 
fact that the Feynman propagator does not solve the homogeneous Dirac equation. 
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Therefore Dencker's theorem on the propagation of singularities can only be used to 
determine the polarization vectors away from the diagonal. Nevertheless, we expect 
that the Feynman propagator for physical states does not have a larger polarization 
set than on Minkowski spacetime. The knowledge of the singularity structure of the 
Feynman propagator should then contribute to a better understanding of interacting 
quantum field theories containing fermions. 
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